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Abstract 

A new space-time model for interacting agents on the financial 
market is presented. It is a combination of the Curie- Weiss model and 
a space-time model introduced by Jarpe Properties of the model 
t "^j H ■ are derived with focus on the critical temperature and magnetization. 

O" 1 It turns out that the Hamiltonian is a sufficient for the temperature 

parameter and thus statstical inference about this parameter can be 
performed. Thus e.g. statements about how far the current financial 
situation is from a financial crisis can be made, and financial trading 
stability be monitored for detection of malicious risk indicating signals. 



g : 1 Introduction 



The foundation of the Ising model became an very important event in the 
modern physics. It was the basic tool explaining critical temperatures for 
which phase transitions occur in physical systems (see Domb et al [1]). It is 
one of the most studied models with wide range of applications in different 
sciences. Weidlich used Ising model to explain the polarization phenomena 
in sociology. It also has been adapted in economics to explain the diffusion 
of technical innovations. New technologies were stated as the result of the 
interaction with neighboring firms. Into general equilibrium economics Ising 
model was introduced by Follmer [2]. There are several sources of influence 
on the price of the stock. One important component is real firm data; another 
one, correlation between amount of buyers and sellers. The action of each 
trader (he is a buyer or seller) was taken as the value of the trader. Kaizoji [5] 
introduced an Ising-type model of speculative activity, which explain bubbles 
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and crashes in stock market. He introduced the market-maker, who adjusts 
the price on the market in dependence of correlation of the buyers and sellers. 

After Curie had discovered the critical temperature, Weiss developed a 
theory of ferromagnetism based on a spin system. It appears by replacement 
of the nearest-neighbor pairs interacting of the Ising model by assumption 
that each spin variable interacts with each other spin variable at any site of 
the lattice with exactly the same strength. 

Some space-time models have been suggested. One model based on the 
Ising model, was suggested by Jarpe |Q. The partition function of the model 
contains two Hamiltonians, one of which describes previous moment of the 
time. This model may be used for describing volatility clustering on the 
market. 

The purpose of this paper is to develop a new space-time model, which is 
also describes the volatility clustering without assumptions about the struc- 
ture of the lattice. Today information about trading on the markets is avail- 
able, e.g. in Internet, to everyone. Since there is less space restrictions for this 
reason it motivates creating a space-time model baced on the Curie- Weiss 
model. Such a model is formally defined in this paper, and some results about 
critical temperature of the market is derived according to the distribution of 
the Hamiltonian in this model. 

In Section |2] the model and methods which are used are described. In 
Section [3] all main results are presented. The implications are discussed in 
Section HI 

2 Model and Methods 

The model which is introduced in this paper is based on two models. The first 
is the Curie- Weiss model, which is a simple modification of the Ising model. 
It allows all agents in the system to interact with each other with a constant 
strength. The second model is the spatio-temporal model of Jarpe [I] which 
possesses both spatial and time dependence. The state of a site in a lattice 
is depending on the states of its nearest neighbors and on the global degree 
of clustering of the previous pattern. 

We took from the Curie- Weiss model the idea of the global interaction and 
from the model of Jarpe the structure of the partition function and obtained 
a new space-time model which is appropriate for describing the volatility 
clustering on the market. 

Let us consider a market which contains N traders symbolically denoted 
by % = 1, 2, . . . , N. In this simple model every trader in a time-period can buy 
a fixed amount of stock or sell the same amount. In the first case the " Trader's 
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decision" is X; = 1, in the second case X; = — 1. Then X = (X 1 , X 2 , ■ ■ ■ , X N ) 
represents the investment attitude of the market. All traders are neighbors. 
That means that each trader knows about the "Trader's decisions" of all 
others traders, so his decision is under influence of the others. 

A configuration of the model is a specification of "Trader's decisions" of 
all traders of the market. With each configuration x — {xi : % — 1, 2, . . . , N} 
a Hamiltonian or interaction energy, 

^ N-l N 
i=l j=i+l 

is connected. We will consider this Hamltonian without investment environ- 
ment. 

Let p(x k ) = P(X = x k ) be the probability of observing the state x k where 
(x 1 , x 2 , . . . , x 2 ) is an enumeration of the distinct states of X. Obviously we 

have that < p(x h ) < 1 and Ylt=iP( xk ) = 1- Further we assume that all 
states x are possible (i.e < p(x k ) < 1). 

Now, wanting to minimize the entropy we have an optimization problem 
of minimizing 

2 jv 

E P (H{x)) = J2p(x k )H(x k ) (1) 

k=l 

with respect to measure P. 

Suppose that the energy of each configuration x^ has been determined. 
The probability, P, that the system has configuration x with energy H if the 
configuration at the privious moment of time is given is: 

p(x\x') = P(X t = x | X t _ x = x') = Z- 1 exp(§H(x)H(x')) 

where Z x > = Y2 X exp(^H(x)H(x')) is a partition function and (3 is the market 
temperature describing the strenght of interaction between the traders. 

2.1 Sufficient statistic 

Proposition 1 The statistic H(X t ) is minimal sufficient for inference about 
the temperature paramter conditional on the previous state, H(X t _ 1 ). 

2.2 The Critical Temperature Of The New Time-Space 
Model 

The behavior of the system in the Curie- Weiss model is described by the 
equation 

m = tanhm (2) 
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where m(x) = -k Yli x i represents magnetization of the configuration. This 
equation allows us to obtain the property of the temperature of the market 
and the critical value of the temperature. 

We will use the method of Hartmann and Weigt [3] to obtain the equation 
of the behavior of new model. 

Theorem 1 The equation of the behavior of the system for the new model 
is 

m(t) = tanh — ^ ^ '— (3) 



Theorem 2 The critical temperature of the new model is (3 = 2. For < 2 
the model {X t : t G Z} is stationary in the space-time sense. For (3 > 2, the 
variable X t is stationary in the spatial sense conditional on X t _i = x t ~\ if 



\m(t-l)\ < v ... 

When N — > oo the partition function has the form 



exp(-Nf{m {t),m(t-l))), 
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Nf"(m (t),m(t-1)) 
We know the function / and the form of the /" is deduced from 
f(m(t),m(t-l)) = l±^ln(i±fW) + l^ln(^fW)-fm 2 (t)m 2 (t-l) 
implying that 

f(m{t),m{t- 1)) | m(t)=0 = ln± 
r(m(*),m(t-l)) | m(t)=0 = l-fm 2 (t-l) 

Therefore 
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N(l - fm 2 (t- 1)) 

Since the Hamiltonian H is sufficient for the temperature parameter, we 
are interested in obtaining the distribution of the Hamiltonian. Testing for 
dependence will make a null hypothisis assuming independence, and thus we 
first consider the distribution of H assuming (3 = 0. We are interested in 
analysis of dependence between and Xj for i,j = 1, . . . , N. 

Theorem 3 Let the variables {Xi : i = 1, . . . , iV} be independent of each 
other and take their values in { — 1, 1} with equal probability |. Then all non- 
identical pairwise products are independent, i.e. XiXj-LXkXi if i =j,k,l or 
j = i,k,l for any dimension N. 
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In this paper, we consider a model where a and b are decisions of the 
traders. For every trader the probability of their decision to 'buy' or 'sell' is 
| unconditional of the other traders. 

Now, let us consider the set of values of the mean field M = ^2 i=1 Xi. 
For this we have the state space 

{-1, -! + #,-! + £,. ..,1} 

meaning that M can be in N + 1 states. Consider now values of M 2 . If iV is 
even, then one value of M is and 

m 2 g{i,(i-|) 2 ,(i-a )2 ,...,o} 

which is to say that M 2 could be in y + 1 different states. If N is odd, then 
is not a value of M and 

M 2 G{l,(l-f) 2 ,(l-|) 2 ,...,^} 

and in that case M 2 could be any of states. 

Regarding the Hamiltonian H(X) = -jfT li<j X i X j = ~^M 2 + \ we 
consequently the statespace of H 

f l-N 1 _ iV/i _ 2_\2 1 _ JV/-1 _ _4 \2 11 
I 2 ' 2 2 V 1 N' ' 2 2 V 1 N> ■>•••■> 2) 

in the case when iV is even. If iV is odd we get 

f l-N 1 _ AT/i _ _2_\2 1 _ iV/i _ J_\2 1 _ _2_ 1 

I 2 ' 2 2V 1 N> ■> 2 2\ L N> ■>■■■■> 2 N 5 

Lemma 1 // P(Xi = 1) = p and {Xi : i — 1, . . . , iV} are independent, then 
«±il G Bin(N,p), where M = i Eii^- 

Theorem 4 TTie distribution of XiXj when {X^} are independent is 

f/VxX -h\ - \ ( ^ Hh=-N/2 
\^ ' >~ ) ~ \ ( (v ^+W^ if h* -N/2 

2.3 Hamiltonian distribution with dependent traders 

Let us now consider the case when the decisions of the traders are not inde- 
pendent. 

Theorem 5 The distribution of Yli<j XiXj when {Xi} are not independent 

Z -^- PI2 {n/2)w ifh=-N/2 



i<j 

is 



S-xX-h\ \ Z e '"^,2)w ifh=-N/2 

h ) ~ l z ' lehm ((vm§ +N)/ 2)w^ ifh^-N/2 
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2.4 Time dependent process 

From now on we consider the space-time process X = {X t : t G Z} where 
X t = {X iyt : % = 1,...,N}. Then we have a corresponding sequence of 
Mean fields, {M t : t e Z} where M t = jjYliLiXi,t, and of Hamiltonians, 
{H t :te Z} where H t = -± ^ £f =m X ht X ht . 



3 Results 

Theorem 6 The sequence of Mean fields, {M t }, is a Markov chain with 
transition probabilities 

r , . .. . ,'(3(1 — NmT){l — Nrrii ' 

P(M t = m t | M t _i = m t _i) 
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Theorem 7 The sequence of Hamiltonians, {H t }, is a Markov chain with 
transition probabilities 



P(H t = h t \H t _ 1 = h t _ 1 ) = 



2 ( (N+ ^m^)/2) z hL e Mih t h t ^) ifh t ^i 

Theorem 8 The conditional expectation of the Hamiltonian is 



and the conditional variance is 

V(H, | = ht-i) = (^.^t± 

3.1 Asymptotics 

Theorem 9 In case with independent trander we have for large N 

Theorem 10 For large N the conditional expectation of the Hamiltonian is 

N(l-2ht-i) 



E(H t | H t _! = ht-!) 

and the varinace is 

V(H t | H t _, = ht-!) 



Ah t -!(N - § (1 - 2h t -!)) 

N\l-2h t ^f 
8/ i 2 _i(^-f(l-2/ it _ 1 )) 2 
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3.2 Stationarity 

The process {X t } is a time-homogeneous Markov chain because 

P(X t+1 = x \X t = x') = P(X t = x | X t _i = x') 
for all states a: and x' and time-points t. 



Theorem 11 The process {X t } is time-reversible and has stationary distri- 
bution 

TT{X) 



where Z x = J2 v ^ x p(n H ( x ) H (v)) 



3.3 Exact calculations 

Example 1 We obtained the exact distribution of the statistic Ylii<j ^i,t^j,t- 
Now let us calculate this distribution in the case with 10 traders. We wrote a 
program in the R language of programming which calculates the probability of 
all possible configurations of the system containing N traders, by determining 
the value of the Hamiltonian for each configuration and build the matrix Vec 
of dimensions 

3 x ("the largest possible value of H" — "the smallest possible value of H") 

where the first line contains all possible values of H , the second line: the 
number of the configurations that leads to the corresponding value in the 
first line, the third: probability of corresponding value from the first line. In 
Figured we can see the distribution of H for N = 10. In the TableUl are the 
exact numerical values of the energy distribution. 



h 


-5 


-3 


3 


13 


27 


45 


Number of occurrances 


252 


420 


240 


90 


20 


2 


Probability 


0.246 


0.41 


0.234 


0.088 


0.0195 


0.00195 



Table 1: The exact distribution values of the energy function when there are 
10 traders. 
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1 1 1 1 

10 20 30 40 

The possible values of H 

Figure 1: Exact distitribution of the energy function when there are 10 
traders. 

3.4 Hypothesis test of independence 

If there is only weak interaction between the traders, then the Hamiltonian 
is more likely to attain smaller values. If the interactions are stronger, then 
larger values of the Hamiltonian is more likely and one talks about mag- 
netization, which could be an indicator of, or even a cause of, bubbles and 
crashes on the market. Therefore methods to state wheter the values of the 
observed Hamiltonian is in some dangerous region is of vital importance to 
decision makers and inderictly to the whole society. 

In order to see if the value of the Hamiltonian deviates from zero to such 
an extent that dangerous development is indicated, the correct distribution 
is needed. If the Hamiltonian approaches the critical value, the system may 
be in a dangerously instable state and a bubble or a crash on the market can 
appear. 

Assume that we make n observations x±,X2, ■ ■ ■ ,x n of the system. Then 
the corresponding Hamiltonian values h±, h,2, ■ ■ ■ , h n are calculated and cat- 
egorized into K classes. If the number of traders are 10 (as in the previous 
example), then the we may have 6 classes corresponding to the 6 possible 
states of the Hamiltonian. But of course we may define these classes in any 
way we choose. 

Strong interactions as opposed to near independence is reflected by the 
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hypotheses 

' H : = 
#1 : /9 > 

To have argument for dependence, i.e. to prove Hi, the statistic 

k=l 



s 



E k ? 



Eh 



may be used, where Ok is the number of observations in klass k, and E k is 
the expected number of observations in class k according to the distribution 
of S under H which is xl-i- Thus the null hypothesis is rejected at level a 
of significance for values of S greater than C which is 1 — a percentile of the 
X 2 distribution with k — 1 degrees of freedom. 

Example 2 Let us consider the data from the Swedish steel market for the 
October 22, 2008. We have the traders and buyers with certain moments of 
trading. To analyze these data we will collect information about trading by 
dividing the time on parts of each about 10 minutes. We than have the ten 
most active traders (AVA, CSB, DBL, ENS, EVL, MSI, NDS, NON, SHB, 
SWB) and twenty intervals of activity. Let us consider the Hamiltonian for 
these traders. The result is in Tabled Thus N = 10 and K = 4. The value 



State 


-5 


-3 


3 


[13,45] 


Number of occurances 


2 


13 


4 


1 


Expected number of occurance 


4.92 


8.2 


4.688 


2.19 



Table 2: The frequencies of observations divided into four classes. 



of the test statistic in this case is 



(2-4.92) 2 . (13-8.2) 2 . (4-4 
4.92 ' 8.2 ' 4.688 



2 , (1-2. 19) 2 
' 2.19 



5.29 < 7,8147 = X3(0-05) 



This means that we can not reject the hypothesis of independence on level 070 
of significance (or any lower level). As a matter of fact the p-value of S in 
this case is 0.15175 so dependence can not be proved on any reasonable level 
of significance. 



4 Discussion 



In this thesis we investigated a new space-time model for interacting agents 
in the financial market. First we reviewed the history of the Ising model 
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and some other Ising-type models. Then the Ising model, Curie- Weiss model 
and some modifications of these models were formally presented. Also we 
considered one way of finding the critical temperature of the market. A 
new space-time model was developed and necessary and sufficient conditions 
for its stationarity were found. The non-linear sensitivity of market global 
properties in terms of temperature parameter changes was investigated. The 
critical temperature for this model was analytically derived. 

The distribution of the Hamiltonian was analyzed using its dependence 
with the magnetization of the market and the exact distribution was cal- 
culated. The conditional expectation and variance of the Hamiltonian were 
found and the stationary distribution was obtained. Then the exact distri- 
bution of the Hamiltonian for 10 traders was calculated, and the expected 
distribution was confirmed. Hypothesis test for independence between agents 
was considered for the Swedish steel market and it showed that there is no 
evident critical situation on the market at the time of this dataset. 

The parameter reflects how strongly traders are influenced by each other 
in the market. It can signalize risk for a crash or a bubble in the market. 
Therefore it is very important that its analysis is accurate in such situations. 
What remains for future work? We could try a lot of real data and compare 
inferential results relying on this model to observable quantities generally 
accepted as a measure of health of the situation on the market. Also we 
can consider a bigger ammount of traders to have a more exact p-value in 
the hypothesis testing. Then we can estimate the amount of interaction for 
our model using e.g. maximum likelihood estimator. Here exists also the 
possibility to develop hypothesis testing based on a time dependent model. 
Also interesting to find out how good this model is to explain volatility 
clustering. 

References 

[1] Domb, C, Green, M.S. and Lebowitz, J.L. (2001) Phase Transitions and 
Critical Phenomena. Academic Press. 

[2] Follmer, H. (1974) Random Economies with Many Interacting Agents. 
Journal of Mathematical Economics, 1(1), 51-62. 

[3] Hartmann, A.K. and Weigt, M. (2005) Phase Transistions in Combina- 
torial Optimization Problems. Wiley- VHC. 

[4] Jarpe, E. (2005) An Ising- Type Model for Spatio- Temporal interactions. 
Markov Processes and Related Fields, 11, 535-552. 



10 



[5] Kaizoji, T. (2000) Speculative Bubbles and Crashes in Stock 
An Interacting- Agent Model of Speculative Activity. arXiv.org. 



11 



